The Dielectric Permittivity of Crystals in the Reduced Hartree-fock Approximation

In a recent article (Cancès, Deleurence and Lewin, Commun. Math. Phys. 281 (2008), pp. 129–177), we have rigorously derived, by means of bulk limit arguments, a new variational model to describe the electronic ground state of insulating or semiconducting crystals in the presence of local defects. In this so-called reduced Hartree-Fock model, the ground state electronic density matrix is decomposed as γ = γ per + Qν,εF , where γ 0 per is the ground state density matrix of the host crystal and Qν,εF the modification of the electronic density matrix generated by a modification ν of the nuclear charge of the host crystal, the Fermi level εF being kept fixed. The purpose of the present article is twofold. First, we study more in details the mathematical properties of the density matrix Qν,εF (which is known to be a self-adjoint Hilbert-Schmidt operator on L(R)). We show in particular that if R3 ν 6= 0, Qν,εF is not trace-class. Moreover, the associated density of charge is not in L(R) if the crystal exhibits anisotropic dielectric properties. These results are obtained by analyzing, for a small defect ν, the linear and nonlinear terms of the resolvent expansion of Qν,εF . Second, we show that, after an appropriate rescaling, the potential generated by the microscopic total charge (nuclear plus electronic contributions) of the crystal in the presence of the defect, converges to a homogenized electrostatic potential solution to a Poisson equation involving the macroscopic dielectric permittivity of the crystal. This provides an alternative (and rigorous) derivation of the Adler-Wiser formula.